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I. ONLINE MATERIAL IN THEORY PAPERS IN NATURE AND SCIENCE. 



The past decade has seen the rise in importance of high impact science journals hke 
Science, Nature and its spectrum of associated journals, Nature Physics, Photonics and 
Materials. Funding agencies increasingly look to measure physicists' performance by the 
articles they have published in these high impact journals. 

The established format for papers in these high impact journals tends to minimize the 
number of mathematical equations, favoring a more conceptual and richly colored figure- 
based representation of key results. This format is ideal for experimental papers, but puts 
theoretical papers that rely on the language of mathematics at a disadvantage. We believe 
that the supporting online materials that accompany Science and Nature articles can provide 
a new format that can help redress this balance. More mathematical theory papers that are 
submitted to these journals can now be written with the main conceptual results in the body 
of the paper, accompanied by key computations and appendix material online. Of course, 
such material can, in the course of time, be groomed for publication in a longer article, but 
in the mean time, this provides a mechanism for key theory papers to be published in high 
impact journals. 

The online material presented here provides the background material for reproducing our 
results. We have given an introduction to symplectic spins, and a full derivation of the 
the application of symplectic N to the new heavy electron superconductors PuCoGa^ and 
NpPd^Al2. We have also included a brief section describing the test-bed application of this 
same method to frustrated magnetism. 

II. SYMPLECTIC SPINS 

A. N-dimensional Symplectic Pauli Matrices 

Symplectic spin operators form a subset of the generators of the SU{N) group. To 
determine their general form, we simply project out the component S of the SU{N) spin 
generators which reverses under time-reversal, i.e the components S for which eS'^e^ — —S. 
For even N, the fundamental SU (N) spin generators can be written 
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Here, all indices range over [±l,±/c], (excluding zero) where N — 2k is even. The general 
symplectic spin operator is obtained by subtracting the time-reversed SU{N) generator 
eT^e^ from T, SP'^ = Tpi - e[rP^Ye^. Putting [e]^ = where a = sgn(a), then 



This traceless matrix satisfies S^'^ = —e[S^'^]'^e'^ , or S*^^ = —aPS^'^ 



(2) 



Since 



—pqS~'^ we can choose a set of ^(A'"-!- 1) independent generators by restricting p + q > 0. 
As in the case of SU{N) matrices, Hermitian generators can be obtained by either sym- 
metrizing, or antisymmetrizing 5"^* on p and q. The resulting matrices form a set of N 
dimensional symplectic Pauli matrices. 



(3) 



where 



pq 



(4) 



^{N + 1) 



V2 

^ 2 i\p\ = \q\), 

normalizes Tr[(7^(7^] = 2Sab in the same way as Pauli matrices. The Dn — ^ 
component "vector" of matrices a^, where [<Tjv]" = cr^ (a = 1,2 .. .D^) plays the role of 
Pauli matrices for SP{N). As an example, consider N — A where the spinor ip and spin-flip 
matrix e take the form 



In this case, there arc 10 symplectic matrices 



10 
-10 
1 
0-10 







(5) 




where a 



(cr^, (7^, cr^) denotes the three possible choices of Pauli matrix. 



(6) 



B. Dot Product 



Here we derive the "dot product" between two symplectic spins. Any even dimensional 
matrix can be divided up into a symplectic and an antisymplectic part M — Ms+Ma, where 



Ms = —eM^e^ and Ma = eMj^e^ . The symplectic part is obtained by projection, M; 



s 



PM, where PMa = removes the antisymplectic component. Now since MA + eMA'e' = 0, 
it follows that 



PM = ^{M -m^e^) (7) 



In components 



1 
2 



ll^sS', - ^^,4]M'\ (8) 



so that 



= li^sS', - e-,e^s]- (9) 
Now, we can always expand Ms = Yla^a<^N terms of symplectic Pauli matrices, and 
with the normalization Tr[a^a^] = 2Sab, ma = |Tr[cT^M], so PM = l^a^r[a%M]a%. 
Expanding both sides gives 

P;^M'^=^-J2KM'^]a:„ (10) 

a 

(where we have temporarily dropped the label on the ajy matrices) or 

pS = ^E</3<^- (11) 

a 

Inserting ([9]), we obtain an explicit expression for the expansion of the dot product between 
symplectic matrices 

J2^^%U{^nU = [StS^, - (12) 

a 

When used to decouple interactions, the first term leads to particle-hole exchange terms, 
while the second term introduces pairing. This same completeness result is also ob- 
tained by brute-force expansion using the explicit spin representation ([2]), which leads to 



C. Abrikosov Pseudo-Fermion representation 

Antisymmetric representations of symplectic spins are obtainsed using Abrikosov pseudo- 
fermions^. An explicit expression for the symplectic spin operator is given by 

S^' = fm^fp = [flf, - Pqf-J-p]. (13) 
4 



where, as before p = sgn(p). (Note the use of the carat over S^'' to dehneate the quan- 
tum operator from the matrix 5''''^. ) The corresponding Hermitian spin operators can be 
obtained by symmetrizing and antisymmetrizing on p and q, as described in ([3]), writing 
= fi{(^N)ai3fi3- Using the dot product relation ([I2]), we can relate these two forms for 
the spin operator via 

= flfp-pqf^-pf^, = S'' (14) 

so that 

SP^ = S- {ctD,,. (15) 

1. SU (2) gauge symmetry 

To examine the properties of these symplectic spins, it is convenient to introduce the 
pairing operators 

a>0 ae[±l,±k] 

vl> = ^ = ! J2 ^f-»fo^ (16) 

a>0 oe[±l,±fe] 

that describe the creation and annihilation of time-reverse pairs of f-electrons. It is also 
convenient to introduce the isospin vector ^ = (\E'i, \E'2, ^3), where ^1 = + \E'), \E'2 = 
-z(^t _ ^) and ^3 = Ea>o(/l/" - f~J~a) = ^/ - f , which satisfy an SU{2) algebra 
[^a, ^b] = '^i^abc^c- The iuversiou of spin under time reversal ensures that the pair creation 
operators creates a spin-singlet of fermions, i.e \E'^ commutes with the spin operator. To 
see this explicitly, note that the commutator with the creation operator 

[S'\ fi] = fliS'^W, (17) 

where we have used an index summation convention over (3 G [±1,±A;]. The commutator 
with the corresponding time-reversed fermion yields the reversed spirip^ 
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so that 



[S'', «/!/! J = fi K - SZ] Pf-, = 0. (18) 

Thus the pair operator [\E'"'', S^'^] = 0. The importance of this relation works both ways: the 
pair is invariant under SP{N) spin rotations generated , while the SP{N) spin generator is 
invariant under the particle- hole rotations generated by \E'^. 

The odd parity of spin operators under time reversal thus ensures that they not only 
commute with the particle number nj, they also commute with the pair operators 

[rif, = = [^t^ = 0. (19) 

These identities imply that the spin is invariant under continuous Boguilubov transforma- 
tions 

fa — 'ufa + v sgn(a)/l^ SU (2) symmetry (20) 

where |up+|fp = l. In this way, for fermionic spins, time reversal symmetry gives rise 
to an SU{2) gauge symmetry of symplectic spins. This symmetry was first discovered for 
spin 1/2 by Affleck et al^'. The above reasoning extends their work, and identifies the SU{2) 
gauge symmetry as gauge symmetry that survives for the (fermionic) generators of SP{N) 
for all N. 

By contrast, had we carried out the same calculation using dipole operators, "P^^ = 
fpfq+pqf-qf-p, which do not invert under time reversal, we would find that the commutator 
of the operator with the time-reversed fermion does not change sign, so that 

[P'', J = fl [Kl+Kl] Pf-P ^ 0. (21) 

This means that the pair operator \E' is not a singlet under the action of the dipole operators. 
It also means that the dipole operators are not invariant under particle-hole transformations. 
Indeed, if the Hamiltonian contains any dipole spin operators, it no longer commutes with 
the singlet pair operator reducing the gauge symmetry back to a U{1) gauge symmetry. 
To fully expose the SU{2) gauge-invariance, it is useful to introduce the Nambu spinor 

/a= I I , /I = (/l,a/_.), (22) 
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By direct expansion, the dot product of these two spinors is the symplectic spin operator: 

fl-f, = flu + Pqf-pf\ = s'' + (23) 

i.e. 

S"' = fl-h- S„ (24) 

Under the SU (2) gauge transformation, fg — > gfg and f^g'' where g = ( ) is an 

SU{2) matrix, so that ■ g'^ g ■ fq = fp ■ fq is exphcitly SU{2) gauge invariant. 



2. Constraint 



To faithfully represent the spin as an irreducible representation, we need to fix the value 
of its Casimir 5*^. If we compute the Casimir using the completeness relation (fT2l) . we obtain 



= J2^fi{cr%umfi{^%ufs) 

= {flfp){flfs)[SasSp^ + eo^^esp]. (25) 

This expression can also be obtained by directly expanding the unconstrained sum 
I X]p q S^'^S'^^. If we normal order the fermion operators in the second term, we obtain 

= nj{N + 2-nf)- J^^aflfljiPf-pfp), 

a,l3 

= nf{N + 2-nf)- (26) 

where Uf = is the number of fermions and we have introduced the pairing terms 

defined in ([H]). then ^' = - z^a), ^'^ = + ^^2) and = ^3 + f , so that 

= ^(f + 2) + 2M/3-(vI/3)'- 4^t? 

= ^(| + 2)-^^ (27) 
since [^'i,^'2] = 2*^3. Ahernatively, 

l(S' + ^f=jij + l), {J = N/A), (28) 
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where, since N is any even number, j is an integer or half-integer. This useful identity 
generalizes the well-known property of conventional spin-1/2 fermions, expressing the fact 
that the sum of spin and charge fluctuations are fixed. In particular, when the isospin 

is zero ^ = 0, the magnitude of the spin is maximized, ^S^ = j{j + 1). We adopt the 
spin maximizing constraint = in all of our calculations. This constraint imposes three 
conditions: 

^sl^) = {nf-N/2m = 0, 

= J]/a/-a|^)=0. (29) 

The first constraint implies that the state is half-filled, with Uf = N/2. The second and 
third terms express the fact that to obtain irreducible representations of the SP{N) group, 
we must project out all singlet pairs from the state These additional constraints become 
particularly important when we come to examine heavy electron superconductivity, for they 
impose the fact that there can be no s-wave pairing amongst the heavy electrons. In a path 
integral approach, we impose the above constraints through the following term in the action 

-fii'^-T)fa, (30) 

where W = (l^i, W2, W^) is a vector boson field that couples to the isospin r of the f-spinor 
and W^^^ ^Wi± iW2. 

Aside: The three component vector boson that imposes the neutrality on the f-spins bears 
close resemblance to the VT-boson in the electro- weak theory of Weinberg and Salam. Indeed, 
the appearance of charged heavy electrons from neutral spins can be closely likened to the 
Higg's effect that occurs in electro-weak theory. If we combine the constraint field on the 
f-electrons with the potential field acting on conduction electrons, we get 



Hc^Hc^ W^-^fjl^ + W% + W^(+V-a/a (30) 




In the symplectic N description of the Kondo effect, the development of a hybridization in 
the primary screening channel corresponds to a Higg's effect, which leads to the combination 
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Z = ^{W^ + e$) becoming massive, forming a high energy plasmon and 7 = — e$) 

forming the new photon that couples equally to / and c conduction electrons, giving the 
f-electrons charge. 



D. Schwinger Boson representation 

Symmetric representations of the SP{N) group, useful for magnetism applications, are 
obtained in a similar way to the fermionic representation, using Schwinger bosons^, S"" = 
bl^{<7'^)a/3bi3, or using ([2]), we can write the more explicit form 

S'' = biS'J^b, = [b% - pqbl^b^,]. (31) 

However, the bosonic constrain is simpler, as we now show. 
By using the dot-product relationship (fT2l) we obtain 



= J2(bUa%Ubp)ib\ia-^),sbs) 

= {bibf3){b\bs)[6asSi3'r + ea-yesis] 

= {b%){blb.) + aP{bib~p){blM 

= (bibabfsbl) - n,] + ^aPibibljib^pbp) + n;,] , (32) 

where Uh = X]q- b^ba is the number of bosons. The in^ in the two terms cancel one-another, 
while, for Schwinger bosons, the pairing terms inside the second term vanish, 6]j,6^q,q; = 0, 
so the final result is 

S' = {bib^bpbl)=nt{nb + N) (33) 

The Casimir of the representation is thus set by fixing the number of bosons. We choose 
the convention 

Ub = NS, (34) 

where upon 

= N^S{S + 1). (35) 

III. HEAVY FERMION SUPERCONDUCTIVITY 

Here we derive a two-channel Kondo lattice model for the NpPd^Al2 and PuCoGa^ 
heavy electron superconductors and construct the mean-field theory for the symplectic large 



limit of our model. We use this to determine the critical temperature of the uniform 
composite pairing instability in the frame of the symplectic large- mean field theory. We 
derive the Andreev reflection off the composite-paired f-electron and show how the crystal- 
field symmetry determines the structure of the gap. Following our discussion on the mean 
field theory, we analyze the fluctuation corrections to the NMR relaxation rate. 

10) = ira^Ti), 




iri^x) 

f Pu^+ : 5/^ f 

FIG. 1: Virtual charge fluctuations of a model Pu^^ (5/^) Kramers doublet into singlet states. 
The addition and removal of an f-electron occur in channels Fi and F2 of different crystal field 
symmetry. 

A. Construction of the model. 

Our two-channel Kondo lattice model for heavy electron superconductivity assumes that 
the ground-state of an isolated magnetic ion is a Kramer's doublet |Fi(t) containing an odd 
number n of f-electrons (Fig. [T]). In PuCoGa^ and NpPd^Al2, the local moments are built 
out of of /—electrons in the 5/ shell. The Pu^~^ ion in PuCoGa^ is a single f-hole in a 
filled j = 5/2 atomic shell, forming a |5/^) Kramer's doublet with n = 5. The situation in 
NpPd^Al2 is more uncertain, the Curie moment extracted from the magnetic susceptibility 
is closest to that of a 5/^ ion with n = 3. 

We assume that the dominant spin fluctuations occur via valence fluctuations into singlet 
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states 

10) ^ ir..) ^ 1^) 

To illustrate the situation, consider PuCoGa^, where |0) = |/^) is a j = 5/2 f-shell. In a 
tetragonal crystal environment, the sixfold degenerate j = 5/2 multiplet of f-electrons splits 
into three Kramers doublets: {T^ , , Fg}. The Kramer's doublet can be written 

iria) = fljO) (37) 

where /p^ creates an {-hole in one of these three crystal field states. To form a low-energy 
singlet, the strong Coulomb interaction between f-electrons forces us to add a second f-hole 
in a different crystal field channel F2. We assume that this state has the form 



|F2®Fi), = ^ J]sgn(a)/^^^/it^ 



JO). (38) 



In practice, there are many other excited states, but these are the most relevant, because 
they generate antiferromagnetic Kondo interactions. In a conventional Anderson model, F2 
and Fi are the same channel. Hund's coupling forces Fi and F2 to be different, and it is this 
physics that introduces new symmetry channels into the charge fluctuations. The simplified 
"atomic" model that describes this system is then 

Hat = ^o|0)(0| + ^i|Fa)(F(T| +E2|0)(0| (39) 

where Ei < Eq,E2. 

In a tetragonal environment, the three Kramer's doublets are determined by 

fL= E (ra||m)/t^^, (a = ±) (40) 

m6[-5/2,5/2] 

where 



re: 


Jr6± 


- 

~ ^±1/2 






f^ 


= cos/343/2 + sin/5/;^5/2 


(41) 




•'rr ± 


= sin/3/^3/2-cos/?/^5/2- 





Here the mixing angle P fine-tunes the spatial anisotropy of the Vf states (see Fig. [2]). Notice 
how the crystal mixes ±5/2 with the =f3/2 states: this is because the tetragonal crystalline 
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T3/2) 



Tf - I ±5/2) 



±1/2) 






FIG. 2: Showing the three crystal field states (F^, , Fg) with a mixing angle (3 — tt/IO, which maximizes the overlap of the F^ witJ] the nearby ligand 
atoms.Tj' overlaps with the eight out-of-plane ligand atoms, whereas Tj overlaps strongly with the four in-plane ligand atoms. Fg is independent of the value of (j), 
and is a/ways aligned along the c~axis, with minimum overlap with nearby ligand atoms. 

environment transfers ±4 units of angular momentum to the electron. A first approximation 
to the crystal field states is obtained by simply setting /? = 0, so that ~ | =i= 3/2) and 
Ff ~ I ±5/2). 

When this atom is immersed into the conduction sea, the f-orbitals hybridize with con- 
duction electrons with the same crystal symmetry. The hybridization Hamiltonian is written 



Hnybr = E[^r+^^+./r^ + Vr-^-Jr^^ + Vr,4,Jre<r + (H.c)] 



(42) 



where ip^^ creates a conduction electron in a Wannier state with crystal symmetry F . The 
matrix elements of this Hamiltonian between the Kramer's doublet and the two excited 

states are 



(43) 



where a = sgn{a). Thus the removal of an electron occurs in a different symmetry channel 
to the addition of an electron. The projected hybridization matrix becomes 



Hhybr = (wL|0)(ra| + aVrM{^^\i^r,-a + H.c 



(44) 



If we now carry out a Schrieffer Wolff transformation that integrates out the virtual charge 
fluctuations into the high-energy singlet states, where the energy of the absorbed, or emitted 
conduction electron is neglected, assuming it lies close to the Fermi energy, then we obtain 

Hk = - (^i|r^>r,.'4,<,(ria| + j2S'4^_.|ri(7')(ria|^r,-.'a), (45) 

(t',(t=±1 
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where 



Ji 



J2 



Eq — El' " E2 — El 
This Hamihonian can be re-written in terms of spin operators as follows 



where we have dropped potential scattering terms and introduced the notation 

S/ = ^ |ria)cr„^(ri/5|. 

a/3 

for the spin of the Kramer's doublet and 



(46) 



(47) 



(48) 



cr^i (0) = ipl^^a^isipnis, o-^^ (0) = ipl^^a^pipr^iB, 



(49) 



for the spin density at the origin in channel Fi and channel 

If we now generalize this derivation to a lattice, the interaction (l47l) develops at each site, 
producing 



■'V 



(50) 



kcr j 

where Sj is the spin operator at site j and c|^^ creates a conduction electron of momentum 
k. We can relate the Wannier states at site j as follows 

ka kcr 

where 

me[-3,3] 

is the form factor of the crystal field state. The Kondo lattice Hamiltonian then takes the 
form 



2k'- 



Ck'/3 • O/e ^ 



kcr kjk'j 

For pedagogical purposes, we work largely with the model in which the matrices $rk 
are taken to be spin- diagonal, giving rise to a simpler form 

H = $^ek4^Ck. + i Yl ^k,k'(cL'^«/3Ck'/3) ■ S/e'^'^'-'')-^^ 

kcr k,k'J 



(53) 
;/>rkl 

(54) 



13 



where 

./kk' = ^101k02k' + ^202k02k'- (55) 

The resuhs obtained using this model are easily generalized to the spin-anisotropic case by 
restoring the spin indices to the form factors. Lastly, we generalize our model from SU{2) 
to symplectic-A^ by replacing the Pauli spin operators (Ta/3 {(yN)ai3, which we write as 

ko- j 

which in its simpler, spin-isotropic manifestation assumes the form 

H = Y. + i E '^k.k'cLck'.^f (J>^('^'-'^^-^^' (57) 

k(T k,k',j 

It is the large N limit of these lattice models that we have solved in our paper. 



B. Decoupling scheme and SU (2) symmetry 

Here we detail our symplectic-A^ decoupling scheme for the Kondo lattice. To derive 
the decoupling procedure, let us first focus on the interaction at a given site, temporarily 
suppressing site indices j. By applying the the dot-product relation ([12]) on the the Kondo 
interaction = J2r ^kt p7j) . we obtain 

Hkv = -^^N ■ S/ = ^iplai'rpfifsiSasS/sy + e^^e^/j) (58) 

which can be rewritten in the form 

Jr 

where the sum over a and (3 runs over the N = 2k spin indices a, (3 G [±1,±A;]. (An 
alternative way to obtain the same result is to write the interaction as Hkf = 'i^^r'SY and 
expand the expression using (fT3i) .) When we cast this normal ordered interaction inside a 
path integral, we can carry out a Hubbard Stratonovich transformation on both terms, as 
follows: 



= -^J2[ i4MifjM ■■ +(^L«/1 J(/^/-/3^r/3)] (59) 



Hkv - 4aiVrfa + Arfl^a) + H.c + iV (^^W^ArAr^ _ ^^^^ 
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This Hamiltonian has a local SU{2) gauge invariance under the transformations 

(61) 





where g = (^^^, ) is an SU{2) matrix. 

The exact solution of the symplectic large limit is provided by the saddle point where 
Vr and Ar acquire constant expectation values. At first sight, it might be thought that 
the mean field erroneously predicts superconductivity under all circumstances! However, 
provided there is only one channel at each site, the effect of the pairing field in fl60l) at each 
site can always be absorbed by a gauge transformation in which 

VU + ^ V|V|2 + |A|2/a. (62) 

We now restore the site label j to all variables. It is convenient to recast the decoupled 
Hamiltonian that results in a Nambu notation, writing 

j,a>0 

where 



HvK = Yl [(V^kV^,./,.) + {fl)Vr,^rja\ + ^Tr[Vf^.Vr,] (63) 



are the Nambu spinors for the f-electron and conduction electron in channel F, while 



V, 



Vr, Ar, 
Ar, -Vr,^ 

describes the Hybridization in channel F at site j. The summation over a is restricted to 
positive values to avoid overcounting. 

We seek uniform mean-field solutions, where Vr, is constant at each site. In this situation, 
its convenient to re-write the Wannier states and f-states in a momentum state basis, 

where Af. is the number of sites and 




Cka 
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Written in momentum space, the mean-field Hamiltonian is then 

k,o>0 \ / 

We can group these terms into a matrix that concisely describes the mean-field theory as 
follows 



(65) 



k,a>0 



where Vk = Vi^^j;. + V202k- This form of the mean-field theory can be elegantly generalized 
to include the effects of spin-orbit coupled crystal fields by restoring the two-dimensional 
matrix structure to the form-factors $rk; 

Vk — Vi$ifc + V2$2k 



C. Mean Field Theory 

To derive the mean-field theory of the uniform composite pair stated, we must diagonalize 
the mean field Hamiltonian 

ft=H''^') (66) 

with Vk = Vi$ik + V2$2k- A simplified treatment of the theory is obtained by assuming that 
'^'rk = 0rkl are spin-diagonal. A more complete treatment using the general matrix form 
factors is given in (IIII Gp . To examine the uniform pairing state we fix the gauge so that 
hybridization in channel 1 is in the particle-hole channel, with V\ = ivi, Ai = 0, i.e Vi = ivi 
while the hybridization in channel 2 is in the Cooper channel, V2 = and V2 = A2T2. The 
eigenvalues uj^ of i^k are determined by 

det(cul - Hk) =00^ ~ 2aku;2 + ^2^0. 
where we have introduced the notation: 

«k = t^k+ + H4 + > 7k = (ekA-t;L)' + 4(^;ikWk)', (67) 

Vlk = Vi(j)ik, f2k = A202k, fk± = ^Ik ± ^2k- (68) 
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The quantity Ck measures the amphtude for singlet Andreev reflection. This quantity is 
unity for the for the simplified spin-diagonal model. When the form factors $rk contain 
off-diagonal components, the above equations still hold, but with the definitions 



Ck = Tr 



i?eTr 



/(40ik0: 



'2k J 



Tr 



Tr 



*k*2k 



(69) 



The eigenvalues of i?k are given by a; = a;k± and cu — — a;k±, where 



Wk± = ^ak±(ak-7k)^^^- 



(70) 



The quantity 



^'lk^'2kCk 



r=i,2 



(71) 



plays the role of the gap in the spectrum. Quasiparticle nodes develop on the heavy fermi 
surface defined by ek = f^k-/'^ directions where Ak = 0. 

The mean field equations are obtained by minimizing the free energy 

= -NT log[2 cosh(/3^k±/2)] + NK J2 ^ 

k,± 

with respect to A and (vr)^ (T = 1, 2), which yields 
1 

K 



1 tanh(a;k±/2T) / Xa^ - ek(ekA - vl_) 
IT 2^ -~ ' 



k± 



2uj 



'k± 



A/«k - 7k 



0, 



2 tanh(u;k±/2T) / ^ , (ck + A)^ + Aiv^x^s^y 
Ik I 2 ± 



k± 



2a; 



'k± 



1 STa? ^ 



2 tanh(cc;k±/2T) 



2c;- 



2± 



V«k - 7k 
ek - A)2 + 4(t;ikSk)^ 



'k± 



k± — "-^^ \ V 7k 

where we have put s^ — l — c^. In the normal phase either vi or V2 is nonzero, corresponding 
to the development of the Kondo effect in the strongest channel. Therefore, there are two 
types of normal phase with two different Fermi surfaces: 



4 

4 

J2' 



(72) 



Ji > J2, V2 — with spectrum 
<^k± 



^ (^k + A ± ^(6k - A)2 + 4<^ . 



(73) 



corresponding to Kondo lattice effect in channel 1, and 
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• J2 > Ji, fi = 0, with dispersion 




(74) 



corresponding to a Kondo lattice effect in channel 2. 

The two normal phases are always unstable with respect to formation of the composite 
paired state at sufficiently low temperature. 

To illustrate the method, we carried out a model calculation, in which the band structure 
of the conduction electrons is derived from the 3D tight binding model: 



and is a chemical potential. Our choice of the form factors is dictated by the corresponding 
crystal structure of the PuCoGas. We take ^^^^ = $r+k for electrons in channel one and 
$2k = "^r^k electrons in channel two. 

As we lower the temperature, the superconducting instability develops in the weaker 
channel. The critical temperature for the composite pairing instability is determined from 
equations fl72p by putting V2 = 0~^. From the third equation with logarithmic accuracy we 
have \og{Txi/Tc) ~ I/J2 which yields 



signaling an enhancement of superconductivity for Ji ~ J2. 

It is instructive to contrast the phase diagrams of the SU{N) and symplectic large 
limits. In the former, there is a single quantum phase transition that separates the heavy 
electron Fermi liquids formed via a Kondo effect about the strongest channel. In the sym- 
plectic treatment, coherence develops between the channels, immersing the two-channel 
quantum critical point beneath a superconducting dome. This is, to our knowledge, the first 
controlled mean-field theory in which the phenemenon of "avoided criticality" gives rise to 
superconductivity. 

D. NMR relaxation rate 

One of the precursor effects of co-operative interference between the two conduction 
channels appears in the NMR relaxation rate just above the transition temperature. The 



ek = — 2t(cos kx + cos ky + cos kz) — fj, 



(75) 




(76) 
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NMR rate is determined by 

n = — — = -— lim (77 

where / is the hyperfine coupUng constant, uo is the NMR frequency and K^_{uj) is the 
Fourier transform of the retarded correlation function of the electron spin densities at the 
nuclear site: 

poo 

K^_{^)^-i ([^+(O,t),^_(O,0)])e-*di (78) 
Jo 

At the mean field level, N ^ oo, the NMR relaxation rate follows a Korringa law. 

Corrections to Korringa relaxation appear in the l/N"^ corrections to the mean field. To 
simplify our discussion, we assume that the Kondo exchange constants are almost degenerate 
Ji ~ J2. In the approach to the superconducting transition, at T > Tc, in principle, we 
need to examine the effects of fluctuations in the hybridization and pairing amplitudes in 
both channel one and two. The anomalous NMR effects we are interested in come from 
the fluctuations in the composite pairs, and as such arc driven by the interference between 
hybridization fluctuations in one channel and pair fluctuations in the other. To simplify 
our discussion we restrict our attention to the corrections induced by the interplay between 
fluctuations in the Kondo hybridization in channel one and pairing fluctuations in channel 
two. The simplified Hamiltonian for our calculation is then H — Hc + Hq + H2 with 

/ X (79) 

^2 = E (01k/.I+q.K,Ck. + 02k^cLAc,4_k,_. + H.C.) 
k,q;CT 

and He describes the conduction electrons. We ignore the fluctuations in the constraint 
fields, which do not couple to the fluctuations in A and V in the lowest orders that we are 
considering. We also neglect the spin-orbit interaction effects by taking the form factors to 
be diagonal in spin space. Our goal is to compute the corrections to the /-spin correlator 
due to the channel interference, i.e. due to the interactions of between heavy electrons and 
fluctuations of slave fields described by H2- The interference corrections to the relaxation 
rate involve the product 0*j^02ic- 

The relaxation rate will be governed by the /-spin correlations: 

Kff{x-T) = -(fV/|(f,r)/^(f,r)/;(0,0)/T(0,0)e"o'''^'^'')„ (80) 
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where {■■■)c denotes the connected Green's function obtained by perturbatively expanding the 
time-ordered exponential in the high-temperature state where the f-electrons and conduction 
electrons are decoupled at the mean-field level. The leading contribution to the temperature 
dependence of the relaxation rate is governed by the diagram on Fig 1 which describes 
the effect of intersite scattering associated with an electron switching from one symmetry 
channel to another as it hops from site to site. To write down an analytic expression for the 
diagram (Fig. 1) we employ the Matsubara correlation functions for the /— and c— electrons 
together with the correlation functions of the slave fields Ky^r.r) = — (TVF(r, r)V"''(0, 0)) 
and Kj\{r,T) = — (TVA(r, r)A"''(0, 0)). Contrary to the situation of only one conduction 
channel, this contribution enhances the screening of the local spins rather then suppressing it. 
As in the case of the composite pair superconductivity, contribution on Fig. 1 originates from 
an interference between the electrons in two conduction channels undergoing a Kondo effect. 
To compute the relaxation rate we take into account for the onset of the Kondo screening 
produces renormalization of the slave boson propagators. In the Matsubara representation 
these propagators are: 



n -1 



J2 



-1 



(81) 



where IlvAi^^) polarization bubbles associated with hybridization fluctuations in 

channel 1 and pair fluctuations in channel 2, which describe the renormalization due to the 
Kondo scattering. The analytic expression for the diagram on Fig.l reads: 

^ ^ ^ Gf{p~kv)Gf{p + q-kv)}Cv{kv) 

ie,iQv,i^A q,P,kv,kA (82) 

X Gc{p + q)Gc{j))lC^{k^)G f{k^ - p - q)Gf{k/^ - p), 
where we employed the four vector notation p = {p, iuo) and included the form factors into 
the deflnition of the conduction electron propagators. The Matsubara frequency summations 
can be performed by employing the spectral function representation for the correlators in 
expression ( l82l) . For example, 

oo 

GUVM= / --^^ (83) 



-oo 



where pf^cip,^) are the corresponding spectral functions. In the high temperature phase, 
where there is no expectation value to the hybridization V or pairing fleld A, the fluctu- 
ation propagators are independent of momentum Kv,/\{p]ifl) = Kv,/^{i^). The resulting 
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FIG. 3: Channel interference contribution to the NMR relaxation rate originating from the inter- 
action of nuclear moments with /-spins in the normal state. Solid lines are the conduction band 
propagators, wiggly lines are the propagators Ky and K/^ of the slave fields for the channel one 
and two correspondingly. The dashed lines are the /-electron propagators. 

expression for the relaxation rate can be compactly written as follows 

oo 

1 1 .duj 



I iy/.(^)/CAH/Cy(-^)^, (84) 



where Wfduj) is proportional to pf{uj)pc{uj). The integral ( !84l) is dominated by the fre- 
quency region near the Fermi surface. Finally, approximating the slave boson functions 
with lCv,/^{oj) ~ Ji^2/ log[(^ ~ i^)/TK^ we obtain the following estimate for the relaxation 
rate 

' ' ' (85) 



TiT iV2log'(r/r;,) + 7r2' 
Our result for the relaxation rate shows an upturn in (TiT)^^ with decrease in temperature, 
in agreement with experimental data of Curro et al.-. 

E. Composite pairing 

Ostensibly, our mean-field theory is that of a two-band BCS superconductor, with hy- 
bridization processes that pair the heavy electrons, and Hamiltonian described by 

«(k)=H^^) (86) 

However, hidden beneath the hood of theory is the underlying gauge invariance that main- 
tains the neutrality of the f-spins. To understand the pairing, we must look not to the 
hybridization pairing terms, which are gauge dependent, but to gauge- invariant variables in 
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the theory. Indeed, it is not possible to say whether the pairing is channel one, or in channel 
two. In the gauge we have chosen, the scattering in channel one is "normal" Vi = ivi and 
pairing takes place in channel two V2 = ^2'r2- But suppose we make the gauge transforma- 
tion flMl) with Qj = —iT2, then 

Vk = if Ik + V2kT2 — > iraVk = f ik^^a - iv2k, 



W = Arg 



ZTaWf-zra) 



(87) 



which transforms the Hamiltonian to one which is now pairing in channel one, and "normal" 
in channel two. The only gauge-invariant statement that we can make, is that supercon- 
ductivity is not a product of one channel or the other, but instead derives from a coherence 
between the two channels. 

To see this, we must combine the gauge- dependent order parameters 



Vi 



V2 




(88) 



into the gauge-invariant composite 



V2Vi + A2^1 l^2Ai-FiA2 
V1A2-I/2A1 V1V2 + A1A2 



Under an SU{2) gauge transformation, Vi gVi, Vl —>■ V\g\ so that V|Vi ^ Vlg^gVi = 
V2V1 is gauge invariant. We shall now show that this matrix is equal to the amplitudes for 
composite pairing and hybridization, as follows: 



J1J2 



il)\{<yN ■ S)V'2 tp\{<^N ■ S)€^l 
ipie'^{(TN ■ S)^2 ipie'^icTN ■ S)e'4)l 



(89) 



In the mean-field theory we have developed, Vi = ivi and V2 = Ar2, so the composite order 
parameter is thus given by 

Ar2 



J1J2 



(iviA2 



(90) 



Thus it is the combination fiA2 that determines the composite pairing that is ultimately 
manifested as resonant Andreev reflection, (see IIII 
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To prove identity (l89ll . we use a path integral approach. Here, it proves useful to employ 
the following matrix representation for the conduction and f-fields at each site j 

fji fj-i ■ 
ft _/-t 

whose columns are made up the Nambu spinors introduced in fl22l) . Using fl23l) . we can write 
iF^Fj)„ = SPi{j) + 6pg, and using ([15]), S^'^ = S ■ (cr^)^,, it follows that F^Fj = S ■ o"^ + 1. 
When we work with a path integral we shall need the normal-ordered version of this result, 

: F^F, :=S-ctI (92) 




In the following derivation we temporarily suspend the site index j of clarity. The no- 
tation introduced in fl9T]) can be used to recast the hybridization terms in the interaction 
Hamiltonian in a more compact form as follows: 



iTaVrfa = 2 XI ^raVr/a = -Tr 

a>0 a 



where Ur = -F\E'p is a two-dimensional matrix operator. In terms of this representation, the 
decoupled interaction Hamiltonian fl63l) at each site assumes the compact form 



f¥^vI 



2 



UrVl (93) 



H 



K 



- E(^ 



+ H.C. + 



N 

2Jr 



Tr[VrV^] 



(94) 



where we have introduced the two dimensional matrices Now to evaluate the expectation 
value of the matrix operator M. = Vi we add a source term to HrK as follows: 



HK[v] = HK + Tr[r]M + {}i.c.)]. 



(95) 



By varying 7], we can read off the matrix Ai, 

6Hr 



6n 



Pa 



(M) 



Now we can combine the source term in (1^3]) with the final trace term in 
Hk[ti\ in the following form 



to rewrite 



H.c. ) + -Tr 
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where 



f 2, 

N 



If we now carry out the Gaussian integral over the Vr, we obtain 



— Tr 
2 



where to leading order in 



J 



J1J2 



-2f] 



N_ 

Ji 



TV 



2 J1J2 



1^1 



So expanding [r/] to leading order in rj, we get 



Hk[v] 



2A^ 



Tr f/rt/^ + 



Ar2 



r]UlUi + H.c 



Differentiating with respect to t] then gives 



J1J2 



(where the final transition from Grassman variables U\Ui to operators requires the intro- 
duction of normal-ordering, denoted by colons). Using the identity (l92l) : F'^F := ajj ■ S, we 
obtain 



ViVi 



J1J2 



J1J2 

'1^ 



where we have taken the transpose of the expression in the last step. Using 



4a 




(96) 



to expand the last term, we obtain 



J1J2 



: [vl/*(o-;vS)vI/^] 



J1J2 



(97) 
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F. Resonant Andreev scattering 

Composite pairing manifests itself through the development of an Andreev reflection 
component to the resonant scattering off magnetic impurities. We can capture this scattering 
in the mean-field theory by integrating out the f-electrons. This leads to a conduction 
electron Green's function of the form 

G(«:)-i = cu-ekT-3-S(«:), (98) 

where k = (k, uj) and 

S(/€) = Vliiu - Ar3)-^Vk (99) 

describes the resonant scattering off the quenched local moments. The hybridization matri- 
ces are written 

Vk = Vil>ik + V2<i>2k = ivi^ik + Al-skTa (100) 

In our earlier discussions, the quantities $rk were assumed to be spin-diagonal. We now 
restore their two-dimensional matrix character, adding a carat to the symbol to denote its 
matrix character. These matrices act identically on particle and hole states^ , commuting 
with the isospin operators r, but they now contain off-diagonal spin-flip terms. 

It is convenient at this stage to examine the off-diagonal structure of the $rk- These two 
dimensional matrices are proportional to unitary matrices, and take the form 

$rk = 0rkWrk 

where is a scalar and ZYpk is a two-dimensional unitary matrix. This matrix defines 
the interconversion between Bloch states and spin-orbit coupled Wannier states. When an 
electron "enters" the Kondo singlet, its spin quantization axis is rotated according to the 
matrix Wpk- When it leaves the ion in the same channel, this rotation process is undone, 
and the net hybridization matrix 

$k$rk = (0rk)' 1 

is spin-diagonal. However, in the presence of composite pairing an incoming electron in 
channel 1 can Andreev scatter from the ion as a hole in channel 2. This leads to a net 
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rotation of the spin quantization axis through an angle about an axis nk that both 
depend on the location on the Fermi surface, as follows 

*2k • *ik = 02k0ik [ck + isk(nk • cr)] 

where Ck = cos(Ck/2), Sk = sin((^k/2). 

Armed with this information, we now continue to examine the resonant scattering off the 
composite-paired Kondo singlet. When we expand the self energy, we obtain a normal and 
Andreev component, given by 



E(k) = EAr(K) + S^(«;). (101) 



where 



2 .,2 
Ik , '^2k 



LO — At3 ijJ + Ats 

and we denote Vik = Vi0ik, ^2k = ^202k- By contrast, the Andreev terms take the form 



(a;2 _ A2) 

Notice how the Andreev scattering contains two terms 



2t'ik^^2k r_^^^^^ ^ a;sk(nk ■ cr)r2] (103) 



a scalar term j^^^\c\fT\ that is finite at the Fermi energy (a; = 0), with gap symmetry 

Ak OC IV$Jk*lk ~ 01k02kCk 



(A^-u;^) 

of the form 



• a "triplet" term —^"^^^ [•Sk(nk • cr)72] which is odd in frequency and vanishes on 
the Fermi surface. 

In practice, the nodes of the pair wavefunction are dominated by the symmetry of the 
function Ck. When an electron Andreev reflects through one hybridization channel into 
the other, it acquires orbital angular momentum. For example, the "up" states of the 
~ I —3/2) and ~ | +5/2) differ by / = 4 units of angular momentum, so the resulting 
gap has the symmetry of an Z = 4 spherical harmonic, or g— wave symmetry. By contrast, 
the up states of the F^ ~ | — 3/2) and Fe ~ | + 1/2) differ by Z = 2 units of angular 
momentum, and the resulting gap has the symmetry of an Z = 2 spherical harmonic, or d- 
wave symmetry, as shown below: 
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r+ ® 





^^^^^ 



FIG. 4: Showing the three possible gap functions: r^" Ty , (gi Fg, and F^ Fg for mixing angle 
/3 = vr/lO. Positive nodes are colored red, while negative nodes are blue. F^" F^ is g-wave(l = 
4), while F^ Fg are both d-wave(l=2). 

G. Dispersion in the presence of strong spin-orbit coupling 

To develop a mean-field theory in the presence of spin-orbit scattering, we need to diago- 
nalize the the conduction electron Green's function. The eigenvalues are determined by the 
condition 

det[ujl-nik)] = 
If we integrate out the f-electrons, this becomes this becomes 

det[cul - n{k)] = {uj^ - X^)Met[g{Ky'] = {u^ - X^)Met[uj - ekTg - 

Now since E(k) oc ^2^^^ ' convenient to factor this term out of the determinant, so that 



det[cul - nik)] = {uj^ - \^)-Met[{iu^ - A2)6?(/t^-^i 



Now 



where 



(cu2 _ \^)g(n)-^ = coiA + D{(T ■ nk)r2) - Bt, + Cn 



(104) 



A 
B 
C 
D 



2 \ 2 2 

e^iuj^ - A^) + Xvl 



(105) 



and vl^ = vli^ ± 



27 



Now if we project the Hamiltonian into states where (nk ■ cr) = ±1, we can replace 
A + D{(T ■ nk)r2 ^ A± Dt2, i.e 

det[ul-n{k)] = [[ _ ;^2) (106) 

The presence of the u"^ — terms in the denominator results from integrating out the f- 
electrons. In actual fact, there are no zeroes of the determinant at = ±A, and the the 
uj'^ — denominators in these expressions act to factor out the false zeros u = ±A in the 
numerator that have been introduced by integrating out the f-electrons. If now expand the 
numerator: 

det[uj{A±DT2) - Bts + Cti] = [lj^ A^ - - - C^] 
= [(^2 _ ^2 _ ^2 )2 _ (2v,^^2^.s^f] - [t^.{uJ' - \') + \v'_Y - [2\v,^V2^c^f . (107) 

Notice that we get the same result for both ±D. Now we know that there is a factor (ci;^ — A^) 
in this expression, so we can write 

det[iu{A ± DT2) - Bts + Cn] = (cu^ - A^) [co^ - 2uj^ak + 7k] (108) 

By a direct expansion of this expression and a comparison of terms with f ll07p , we are able 
to confirm that this factorization works, with 

ftk = vl^ + ^(A^ + e^) 

7k = (ekA - vl_y + (2t;ikt;2kCk)' (109) 

Thus 

det[cul -7^(k)] = [cu^-2cjV + 7k]^ (110) 

The surviving yet crucial effect of the spin-flip scattering is entirely contained in the Ck factor 
in 7k. The Boguilubov quasiparticles in the composite paired state preserve their Kramer's 
degeneracy, with dispersion given by 

as described in section (IIII Cp . 
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H. Crystal Fields determine the gap symmetry. 



Here we calculate the form factors for the two-channel Kondo model in a tetragonal 
crystal field environment. In a tetragonal crystal field environment, the Kramer's doublets 
are given by |r±) = \jm){jm\r±) , [j = 5/2), where from (1411) 



;iii) 



re: 




= l±l/2) 






= cos/3|t3/2) + sin^|±5/2) 






= smp\TS/2) -cos/3|±5/2) 



The matrices representing these crystal field states are then 



a : 



1 0\ + 
^0 1 

sm{P) cos(/5) 
cos{P) sm{p) 
cosiP) sin(/?) 
sin(/5) cos(/3) 

while the overlap of the Bloch states with the spin-orbit coupled Wannier states \jm) {j 
5/2, m G [-j, j]) is given by 



{TQa\jm) 
{T7+a\jm) 
{T7-a\jm) 



112) 



— TfiCT 

Om|ka) = aW^^F^ .k) 



(113) 



cr 



where we have used the Clebsch Gordon coefficient {jm\lm — |cr, | |(t) = criy ^ , 
±1). We obtain the form factors of the Wannier states by multiplying matrices from flll2p 
with the matrix elements (11121) : 



[$rk]aa = ^{Ta\jm){jm\ka) 

m 

The form factors are then given by 



(114) 



VEcY^^i^) + sYiik) ~V6sYi{k) - V2cY^^{k) 
VQsY^^ (k) + V2cY^^ (k) - VEcYi (k) - sY^2 (k) 



;il5) 
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1 



\/hsY^^ (k) - cYi (k) ^fQcYi (k) - ^sY^^ (k) 



(116) 



1 

7f 



v^F\(k) -v^yo'(k) 

where we use the shorthand c = cos(/3), s = sin(/5) to denote the cosine and sine of the 
mixing angle. Each of these functions is time-reversal invariant, namely they satisfy 



rk 



There are basically two symmetry classes of superconductor that are possible in our model, 
one formed from Fg, the other formed from F^ ® F^ (Fig. |4]). The latter is argued 
to be preferably, because it is these two states that have the maximum overlap with nearby 
ligand atoms. The form factors 0rk are determined from 



^ik)^ = ^Tr[$t,<|.ik] 



;ii7) 



For general j3 , these are form-factors with point nodes along the c-axis, but no lines of 
nodes. The gap function is however determined from the symmetry of 



01k02kCk 



4 



Tr[$^k$ik + H.c. 



In our model, we have chosen Fi = F^ corresponding to the out-of-plane ligand atoms and 
F2 = F7 corresponding to the in-plane ligand atoms. In this case, when we expand this 
function in detail, we find it has the form 

0ik02kCk = cos(2/5)Agi(k) - sin(2/?)Ag2(k) 

where A^i and Ag2 are g-wave gap functions of the form 



A.i(k) 



75 
levr 



cos (40) sin^[6'] 



and 



A,2(k) 



167r 



sin2(^)(3 + cos(20)). 



For small /3, the gap is dominated by Ac,i(k) 
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IV. FRUSTRATED MAGNETISM 



As a test of symplectic-A^, we would like to return to the origins of SP{N), initially 
developed to treat antiferromagnetism on frustrated lattices^. SP{N) describes antiferro- 
magnetism, or the formation of valence bonds well, but cannot handle ferromagnetism, the 
fluctuations of those bonds. This weakness can already be seen in the simplest model in 
frustrated magnetism, the Ji — J2 Heisenberg model on a square lattice^ 

H = Ji S'x • 'S'x+^ + J2 Sy_ ■ S:x.+fi', (118) 

where Ji and J2 are the first and next nearest neigbor couplings. 

For large J1/J2, the ground state is the typical Neel state, where diagonal spins are 
ferromagnetically aligned(see Fig [5] inset). As J2 increases, these diagonal bonds become 
increasingly frustrated, destabilizing the long range order, which requires higher and higher 
spin, Sc for magnetic order. 

At large J2/J1, this model describes two interpenetrating Neel lattices that are classi- 
cally decoupled. When fluctuations are included, a biquadratic interaction locks the two 
sub lattices together in a colhnear configuration^(see Fig [5]). The nearest neighbor bond, Ji, 
initially stabilizes the colhnear state, but, as the frustration increases, the long range order 
is eventually destabilized. 

In quantum magnetism, one generally uses a Schwinger boson spin representation. The 
SP{N) approach introduced by Read and Sachdev^^ decouples the antiferromagnetic inter- 
action in the following manner 

Hi = -^Bl,B,2 SP{N). 

where = creates a valence bond (2,1) between sites one and two. The 

SP{N) approach captures the competition between first and next nearest neighbor links for 
valence bonds^, but misses the frustrating effect of the ferromagnetic bonds, resulting in 
overstabilization of the colhnear state^^. If we decompose Hj in terms of spin generators, 
we find it contains a mix of "spins" and "dipoles" 

Hj = ^{S,-S2-Vi-V2), SPjN) . (119) 

This inadvertent inclusion of dipoles with a negative, i.e ferromagnetic sign, tends to can- 
celling out the frustrating effect of ferromagnetic bonds. For instance, in the Ji — J2 model, 
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FIG. 5: The Ji — J2 Heisenberg model: We compare the critical spin Sc = (^)^) below which there 
is no long range order in the ground state, calculated within SP{N) (bold red line), symplectic- 
N (blue and green lines), and spin wave theoryS (thin black line). Ji and J2 are nearest and 
next nearest neighbor antiferromagnetic bonds, as shown in the figure. For small J2/J1, the spins 
configurations are staggered, while for large J2/J1, the ground state breaks lattice symmetry to 
develop collinear order as shown in the figure. SP{N) (bold red line) tends to overstabilize the 
long range ordered phases, most dramatically on the one sublattice side, where the critical spin 
is independent of the strength J2 of the frustrating diagonal bonds^'^. Symplectic- restores the 
frustration-induced fluctuations by treating both ferromagnetic and antiferromagnetic bonds, on 
equal footing, which corrects this overstabilization. 

the critical spin for developing long-range antiferromagnetism is artificially independent of 
the frustration in the the SP{N) mean-field theory^ (Fig- El)- 

In the symplectic- approach, we decouple the interaction exclusively in terms of the 
generators of SP{N) S^p = bl^bjj — a/36^^6Q,. Using the explicit form of the symplectic spins. 
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Hj-^Si-S,--^ 



we find the Heisenberg interaction decouples into two terms 

bI^B2i - AI^A2i , Symplectic A^. (120) 

where A21 = Xlo- ^L^io-- The second term in this interaction describes ferromagnetic correla- 
tions, which were absent in the original applications. The operators ^421 and A\i "resonate" 
a valence bond linked to a third site, between sites one and two: (1,3) ^ (2,3). When 
the bond resonates between sites 1 and 2, the amplitude for singlets to form between the 
two sites is reduced, giving rise to a ferromagnetic correlation between sites 1 and 2. The 
exclusion of dipole spins requires that we treat these two terms in equal measure. 

When we carry out a Hubbard Stratonovich factorization of the Heisenberg interaction 
(11201) . it separates into two amplitudes h and A describing bond resonance and condensation 
respectively 

A 
A h 



JSi- 82= {h\^,ah2 



\+^m'-m- (121) 



This kind of decoupling scheme was first proposed by Ceccatto, Gazza and Trumper— for 
SU (2) spins, where it is one of many alternative decoupling procedures. In symplectic- A^, it 
is the unique form preserving the time reversal parity of the spins. Now we would like to see 
if and when the ferromagnetic h bonds develop and what effect they have on the physics. 
To do this we examine the action, 

NS[h, A, h,\]= / drYSZ. iAdr - \)h. + \NS] + ^ J,,S, ■ (122) 

i a (ij) 

with JijSi ■ Sj given above. We assume all bond fields are uniform and static, depending 
only on i — j. The constraint, rib = NS is enforced by the Lagrange multiplier Aj. As the 
action is quadratic in the Schwinger bosons, they can easily be integrated out to find the 
mean field Free energy: 

m = Jf^ log|2 sinh ^1 + E ^''^<' - - A(2S + 1) (123) 

k 

where («,j) is a pair of sites with nonzero Jj^, M is the number of sites, uj^ = 



a/|A — /ifcP — |AfcP, and /ifc and A^ are the Fourier transforms of hij and A^j. The Neel 
state is described by A's along the nearest neighbor bonds, and induced hd along the diag- 
onal bonds for finite J2- 



^f"' = V + 4/^rfC,c,)2 - 4A2(s, + syf), (124) 
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where q = cos kia and si = sinfc/a. For large J2, the classical decoupled state consists of 
along the diagonal bonds, where the magnitude of is the same on both sublattices, but 
the phase between the two is free. When fluctuations are introduced, the lattice symmetry is 
broken by the collinear order; we choose the phase so the collinear state is antiferromagnetic 
along y{Ay), which induces ferromagnetism along x{hx), giving the dispersion. 



The values of h, A, and A are chosen to minimize the Free energy using the mean field 
equations: dF/dX = 0, dF/dh = and dF/dA = 0, for all h's and A's in the particular 
problem. SP{N) proceeds similarly except that the /I's are fixed to be zero. 

We are interested in the zero temperature case, where in two dimensions, the Schwinger 
bosons can condense when the gap in the spectrum vanishes, signaling the onset of long 
range magnetic order. As we are primarily interested in the spin above which the system 
orders magnetically, Sc, we only consider the unpopulated condensate. 

The minimum gap in the spectrum is fixed at (±7r/2, ±7r/2) and (0, ±7r/2) for the Neel and 
Ising states, respectively. Setting that gap to zero gives us an algebraic relation between the 
parameters supplementing the mean field equations. After the parameters are determined 
from dF/dh and dF/dA, dF/dX = can be used to find the critical spin. 



Results from these calculations for both symplectic-A^ and SP{N) are shown in Fig ??c, 
along with a comparision to spin wave theory. SP{N) drastically overestimates the stability 
of the ordered states, which is corrected by the ferromagnetic bonds included in symplectic- 
A^. For both large A^ theories, the regions of Neel and Ising order overlap, indicating a first 
order transition, while spin wave theory predicts a second order transition for all S, with 
an intervening region of spin liquid. However, when higher order corrections are taken into 
account, modified spin wave theory gives exactly the results of symplectic-A"— . 

We can draw an analogy between frustrated magnetism and heavy fermion supercon- 
ductivity, in which previous large A^ techniques were able to treat only one of two possible 
phenomena, ferromagnetism and antiferromagnetism in the case of frustrated magnetism, 
Fermi liquid physics and superconductivity in heavy fermion superconductors. In both sit- 
uations, symplectic-A^ enables simultaneous and equivalent treatment of both phenomena 




(125) 




(126) 
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and significantly improves upon the previous results. 
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